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Abstract 

The consequences of site-directed mutagenesis experiments are often anticipated by empirical rules regarding the 
expected effects of a given amino acid substitution. Here, we examine the effects of "conservative" and "noncon- 
servativc" substitutions on the X-ray crystal structures of human recombinant FKBP12 mutants in complex with 
the immunosuppressant drug FK506 (tacrolimus). R42K and R42I mutant complexes show 110-fold and 180-fold 
decreased calcineurin (CN) inhibition, respectively, versus the native complex, yet retain full peptidyl prolyl isom- 
crase (PPIase) activity, FK506 binding, and FK506-mediated PPIase inhibition. Interestingly, the structure of the 
R421 mutant complex is better conserved than that of the R42K mutant complex when compared to the native 
complex structure, within both the FKBP12 protein and FK506 ligand regions of the complexes, and with respect 
to temperature factors and RMS coordinate differences. This is due to compensatory interactions mediated by 
two newly ordered water molecules in the R42I complex structure, molecules that act as surrogates for the miss- 
ing arginine guanidino nitrogens of R42. The absence of such surrogate solvent interactions in the R42K complex 
leads to some disorder in the so-called "40s loop" that encompasses the substituent. One rationalization proposed 
for the observed loss in CN inhibition in these R42 mutant complexes invokes indirect effects leading to a misori- 
entation of FKBP12 and FK506 structural elements that normally interact with calcineurin. Our results with the 
structure of the R421 complex in particular suggest that the observed loss of CN inhibition might also be explained 
by the loss of a specific R42-mediated interaction with CN that cannot be mimicked effectively by the solvent mol- 
ecules that otherwise stabilize the conformation of the 40s loop in that structure. 

Keywords: calcineurin; immunophilins; site-directed mutagenesis; structure-based drug design; X-ray 
crystallography 



FK506 (United States Adopted Names Council of the American 
Medical Association, Chicago, Illinois [US AN], tacrolimus) is 
a natural product screening lead (Kino et al., 1987) now ap- 
proved for therapeutic use as an immunosuppressant in Japan, 
the USA, Germany, and other countries. FK506, in complex 
with its 12-kDa M r binding protein FKBP12, exerts its immu- 
nosuppressive effects through the inhibition of calcineurin (CN), 
an intracellular Ca +2 -caImodulin-dependent phosphatase (Klee 
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& Cohen, 1988). CN inhibition, in turn, interrupts the induc- 
tion of IL-2 and other T-cell activation events (Friedman & 
Weissman, 1991; Liu et al., 1991). A homologous natural prod- 
uct, rapamycin (AY-22,989; USAN, sirolimus), which was ini- 
tially discovered as an antifungal agent (Sehgal et al., 1975), can 
antagonize the CN inhibitory activity of FK506 (Bierer et aL, 
1990a; Dumont et al., 1990b), even though it is itself an immu- 
nosuppressant by a different mechanism (Bierer et al., 1990a; 
Dumont et al., 1990a). These differences in CN inhibitory ac- 
tivity between the agonist FK506 and the antagonist rapamycin 
in their complexes with FKBP12 (Table 1), were first explained 
within the framework of an elegant model (Schreiber, 1991) that 
focused attention on the corresponding differences in chemical 
structure between the two ligands in their so-called "effector do- 
main" region (see, e.g., Fig. 1). X-ray crystallographic studies 
have provided support for this model (van Duyne et al., 1991a, 
1991b, 1993; Becker etal., 1993; Rotonda et al., 1993; Connelly 
et al., 1994; Wilson et al., 1995), by showing that the effector 
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Table 1 . Biochemical properties of the native 
and mutant complexes* 





Catcineurin 


Reduction 


PK506 


Specific 






vs. native 




PPIase activity 


Mutant 


(nM) 


complex 


(nM) 


(s- 1 jiM-') 


Nativc b 


5.5(1.8) 




0.6 (0.2) 


4.3 (0.4) 


R42K b 


590 (200) 


107.3 


0.6 (0.2) 


3.8 (0.3) 


R42I b 


970 (150) 


176.4 


0.1 (0.1) 


2.5 (0.3) 


R42Q C 


325 (150) 


59.1 


4.3 (2.0) 


3.0 (0.3) 


Native 11 


7.9 (3.0) 




0.4 (0.2) 


2.2 (0.2) 


R42Q d 


850 (250) 


107.6 


1 .7 (0.6) 


1.3 (0.3) 


R42A d 


280 (80) 


35.4 


0.2 (0.1) 


1.1 (0.2) 


Chimera dc 


19 (2) 


2.4 


0.4 (0.2) 


0.57 (0.05) 



51 Summary of published biochemical data for native and mutant 
FKBP12 proteins and their complexes with FK506and calcineurin. The 
inhibition constant for calcineurin by native and mutant FKBP12 com- 
plexes with FK506 is given, along with the FKBP12 inhibition constants 
versus FK506 and the PPIase specific activity of the various FKBP12s 
versus a synthetic substrate. 

h Aldape et a». (1992). 

c Futcr et al. (1995). 

d Yang et al. (1993). 

tf Substitute FKBP12 residues 40-44 (-RDRNK-) with the corre- 
sponding residues (-LPQNQ-) from FKBPI3. 



domains of FK506 and rapamycin do indeed protrude from the 
surface of their respective protein-ligand complexes (Fig. 2A) 
with distinct conformations that might be compatible with CN 
binding and inhibition for the one ligand, but not for the other. 
In turn, those chemical structure elements shared by the two 
gands (Fig. 1) were shown in those studies to constitute FKBPI2 
"binding domains," allowing a rationalization of the reciprocal 
antagonism between the two ligands in terms of their competi- 
tion for a common FKBP12 binding site. 

The effector domain model has retained broad acceptance as 
a first approximation to the complicated problem of immuno- 
suppressive drug design, in part because of its consistency with 
the observed loss of CN inhibitory activity that follows even mi- 
nor variation in the chemical structure of FK506 (Goulet et al., 
1994). By limiting the role of the FKBP12 protein to that of a 
presenter of ligand functionality to CN (Schreiber, 1991; Rosen 
& Schreiber, 1992; Schreiber et al., 1993), the model reduces the 
scope of the drug design problem to one of simple mimicry of 
the conformation of the FK506 effector domain that protrudes 
from the surface of the native FKBP12-FK506 complex. Unfor- 
tunately, these efforts have yet to produce a linear or macro- 
cyclic drug lead — let alone a clinical candidate— that exceeds the 
potency of FK506 (Itoh et al., 1995); ligands predicted on the 
basis of this model have all turned out to be antagonists of FK506 
(see, e.g., Bierer el al., 1990b; Somers et al., 1991; Armistead 
et al., 1995). 

Experimental evidence for a more complicated interaction be- 
tween CN and the FKBPI2-FK506 complex first emerged from 
a systematic examination of the biochemical properties of site- 
directed mutants of charged residues on the surface of FKBP12 
(Aldape et al., 1992). The critical involvement of the "40s loop" 
and "80s loop" regions of the protein (as defined in Fig. 2B) was 
established for this interaction by these and subsequent studies 
(Yang eta!., 1993; Futeretal., 1995), leading to a generalization 




OCH 3 
FK-506 




Rapamycin 

Hg. I. A: Chemical structure of FK506 (USAN, tacrolimus). The ef- 
fector domain of FK506 (Schreiber, 1991) corresponds to those portions 
of the ligand (C18-C23 and substituents in the macrocycle, and C26- 
C34 in the cyclohexyl ring) that have been shown crystallographically 
to protrude from the surface of its complex with FKBP12; structural 
elements in common between rapamycin and FK506 have been shown 
crystallographically to bind the PPIase active site of FKBP12 in the same 
manner (van Duyne et al., 1991a, 1991b, 1993; Becker et al., 1993; 
Rotonda et al., 1993; Armistead et al., 1995; Itoh et al, 1995; 
Wilson et al., 1995). B: Chemical structure of rapamycin (AY-22989, 
USAN, sirolimus). The rapamycin effector domain corresponds to at- 
oms C15-C29 in the macrocycle and C36-C42 in the cyclohexyl ring. 



of the effector domain model in the direction of a composite 
"effector surface" of both protein and ligand structural elements. 
Elsewhere, we have explored the structural consequences that fol- 
low substitutions in the 80s loop region of FKBP12 (Itoh et al., 
1995) and have identified composite features on the FKBP12- 
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Fig. 2. A: Comparison of the struc- 
tures of the FK BP 12 complexes with 
FK506 (in red) and with rapamycin 
(in white). A dot representation of 
the surface of native FKBP12 (Wil- 
son et ah, 1995) is shown in blue and 
is representative of the surface of 
the FKBP 12 complex structures with 
FK506 and rapamycin. B: Conforma- 
tion of the backbone Ca of FKBP 12 
in its complex with FK506 is shown 
in blue. The 40s loop and 80s loop re- 
gions of the FKBP 12 protein are spe- 
cifically identified in yellow and red, 
respectively. The 40s loop is made up 
of residues 40-44 in FKBP 1 2, which 
form a bulge on the third 0-strand of 
the protein, as defined by van Duyne 
et al. (1991a). The 80s loop includes 
residues 84-91 on the edge of the 
Fifth 0-strand of the structure. 



FK506 effector surface that might be candidate CN recognition 
and binding elements (Wilson et al., 1995). In this paper we focus 
on structural elements in the 40s loop region of FKBP 1 2 and, in 
particular, on substitutions at residue 42 of the protein, which 
demonstrate profound (Table I) but complicated effects on the 
CN inhibitory potency of the corresponding mutant complexes 
with FK506. These mutant data have led to proposed mechanisms 
of action, which are distinctly different in their character and con- 
sequences, that need to be resolved (Clardy, 1995). 



Results 

Crystallization, data collection, and refinement statistics for the 
native and mutant complex structures reported here are given 
in Table 2, along with RMS differences in conformation and mo- 
bility versus the native complex structure; biochemical data are 
summarized from the existing literature (Aldape etal., 1992; 
Yanget al., 1993; Futeret al., 1995) in Table I. In all the com- 
plexes studied, the FKBP12 fold (Fig. 2B) that was seen in the 
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Table 2. Summary of the wild- type and mutant 
complex structure analyses* 





Wild type 


R42K 


R421 


Arcu detector used 


Siemens 


Rigaku 


Rigaku 


P4 2 2,2 unit cell; a, c (A) 


58.39, 55.76 


58.31, 55.93 


58.25, 55.98 


Resolution (A) 


6.0-1.5 


6.0-1.5 


6.0-1.6 


No. observations 


76,344 


43.803 


44,188 


<tt> Reflection (/> 2a) 


81.4 


88.1 


89.7 


tf-mcrge (%) 


3.99 


5.88 


3.15 


tf-factor (%) 


16.6 


18.7 


16.8 


No. water molecules 


85 


83 


87 


RMS bond length error (A) 


0.016 


0.016 


0.018 


RMS bond angle error (deg) 


2.74 


2.85 


2.93 


Avg. FK506 factors (A 2 ) 


12.3 


15.3 


11.0 


FK506 RMS diTf vs. wi (A) 




0.139 


0.123 


FKDIM2 RMS diff vs. wt (A) 




0.146 


0.147 



u Native and mutant FK BP 1 2 complexes all share the native crystal form first 
reported by van Duync el at. (1991a). 



native complex structure (van Duyne et ah, 1991a) is strongly 
conserved, consistent with the observed retention of PPlase and 
FK506 binding activity (Table I). In this study, a considerable 
effort was made to crystallize all the complexes reported in a 
common crystal form (Table 2), in order to facilitate a direct 
comparison between structures. This crystal form turned out to 
be that of the native FKBP12-FK506 complex (van Duyne et al., 
1991a), as a consequence of seeding mutant complex crystalli- 
zation experiments with microcrystals of the native complex and 
subsequently using crystals from those solutions as macroseeds 
leading to data quality mutant complex crystals. 

In the native complex structure (red coordinates in Fig. 3A, 
Kinemage I), the two guanidino nitrogens of R42 are seen to sta- 
bilize the "40s loop" of FKBP12 through their participation in 
a bridging network of noncovalent interactions between residues 
D37 and K44. In the R42I mutant complex structure (yellow co- 
ordinates in Fig. 3 A, Kinemage 1), two tightly bound water mol- 
ecules (yellow in Fig. 3A, Kinemage 1) substitute for the missing 




Fig. 3. A: Comparison of the structures of native (red) and R421 mutant (yellow) FKBP12 complexes with FK506. Refined co- 
ordinates (Table 2) of both complex structures arc superimposed in the figure, together with the 2|F„| - \F C \ electron density 
(in blue) of the R42I mutant complex structure, contoured at I a above background. Part of the FK506 binding domain is shown, 
as well as residues surrounding R42 and H87 in the protein, Water molecules bridging residues D37 and K44 are shown in green 
and yellow, corresponding to the R421 mutant complex structure. Dashed lines indicate bond distances between the water mol- 
ecules. Our observations are inconsistent with the suggestion of Yang et al. (1993) that R42 mutants exert their effects indirectly, 
by reorienting nearby regions of the complex structure. Nor does the conformation of the 40s loop change as drastically as that 
of FKBP13 (Schultzet al., 1994) as a consequence of these substitutions. B: Dot-surface representation of the R421 mutant com- 
plex structure in the vicinity of FK506. Two water molecules (shown in green) in the R421 mutant complex structure fit readily 
into the gap created by the R42I substitution and act as surrogates for the missing guanidino nitrogens of R42 in bridging resi- 
dues D37-K44. This interaction helps preserve the native conformation in the mutant complex and maintains PPlase and FK506 
binding activity. Nonetheless, CN inhibition is drastically reduced (Table 1), suggesting a specific protein-protein interaction 
in the native complex that the water surrogates would be unable to mimic. C: Structure of the R42K mutant complex (in green) 
compared to that of the native complex (in red). As above, the 2|F 0 | - \F?\ electron density corresponding to the R42K mu- 
tant complex, contoured at lo above background, is presented in blue. {Continues on facing page.) 
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Fig. 3, Continued. 
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arginine guanidino nitrogens of R42, and are well accommodated 
(Fig. 3B, Kinemage 1) in the gap created by the smaller isoleucine 
substituent; a third water molecule (green in Fig. 3A, Kinemage 1) 
is apparently common to all the FKBP 12 complex structures. As 
guanidino nitrogen surrogates, the two water molecules help to 
preserve the native complex 40s loop conformation in the R42I 
mutant complex, with an RMS deviation of 0.147 A between the 
structures (Table 2). This interaction resembles one seen in the 
structure of a T157G mutant of T4 lysozyme (Alber et al., 1987; 
Matthews, 1993), where an ordered water molecule, acting as 
a surrogate for the missing T157 hydroxyl, preserves the pattern 
of stabilizing hydrogen bonding interactions seen in the native 
protein. It is interesting that, in spite of the high degree of con- 
servation seen in the structure of the R42I mutant complex, CN 
inhibition is nonetheless reduced by -180-fold (Table 1). 

In the structure of the FK506 complex with the more conser- 
vative R42K mutant protein, the single e-amino nitrogen of the 
lysine side chain is unable to substitute for both of the guani- 
dino nitrogens of arginine (red in Fig. 3C, Kinemage 1); nor does 
the substitution leave enough space for an additional water mol- 
ecule to insert itself as a surrogate, as shown in Figure 3B 
(Kinemage 1) for the R42I mutant complex. The resulting de- 
stabilization of the 40s loop is reflected in the higher tempera- 
ture factors observed in that region of the protein (Fig. 4), even 
though the conformation of mutant protein complex still closely 
resembles that of the native (Fig. 3C, Kinemage 1), and FK506 
binding and PPIase activity are preserved (Table 1). 



Discussion 

Two working models have emerged to rationalize the profound 
though complicated effects on CN inhibition that accompany 
substitutions in and around residue 42 of FKBP12. The simpler 
of these is made evident in the R42 single-site mutant complexes 
that were first characterized biochemically by Aldape et al. (1 992), 
whose crystal structures are described here. In that model, the 
observed loss of CN-inhibitory activity can be immediately ex- 





-.Wild typ. 




— R«I 




— R43K 




I I I I I I 1 I 1 I I I I I I I I I I I I 

FKBP12 Residue Number 

Fig. 4. Average X-ray temperature factors for the main-chain heavy 
atoms of FKBP 12 for the native and the two mutant complex structures. 
Residues showing high temperature factors correspond to loop regions 
in the structure of FKBP 12. These have been shown previously to be re- 
gions of local flexibility within the protein that were identified through 
a crystal structural analysis of 19 FKBP12-ligand complexes, each in 
a different crystal-packing arrangement (Wilson et al., 1995). 



plained in terms of a direct and localized perturbation, by the 
substituted residue, of the effector surface presented to CN by 
the corresponding FKBP12-FK506 complex. The Merck group 
(Becker et al., 1993; Rotonda et al., 1993) arrives at a similar 
conclusion in their analysis and comparison of the structures of 
human and yeast FKBP12-FK506 complexes and of the human 
FKBP12 complex with L-685,818, an l8-hydroxy,21-ethyl an- 
alog of FK506. 

Yang et al. (1993), however, have suggested that substitutions 
at R42 exert their influence on CN inhibition indirectly, through 
a generalized conformational misorientation of nearby elements 
of the FKBPI2-FK506 effector surface. This model is inferred 
from the curious pattern of CN-inhibitory activity evidenced in 
the FKBP12/I3 chimeras studied by these workers. Substitution 
in FKBP 1 2 of the corresponding 40s loop sequence from FKBP 1 3 
(i.e., replacing the sequence RDRNK with LPQNQ) leads to 
only a modest loss of CN-inhibitory potency (by -2-fold to 
19 nM). In turn, the single site R42Q mutant complex is severely 
compromised (by -100-fold, to 850 nM), even though the R42Q 
substitution is incorporated in the FKBP 12/ 13 chimera. From 
these observation, Yang et al. (1993) concluded that the effects 
of an R42 substitution would have to be strongly contextual. In 
other words, R42 and Q42 would each be appropriate to the 40s 
loop of FKBP 12 and FKBP 13 respectively, with only a modest 
loss of activity for the latter in the chimeric FKBP 12/ 1 3 com- 
plex. An incompatible substitution, such as that of R42Q into 
FKBP 12, would then lead to a significant and generalized dis- 
ruption of the effector surface, an event that would be reflected 
in the much lower CN-inhibitory activity seen in the single- 
mutant complexes. Clardy (1995) has noted that the 40s loop in 
the structure of the native FKBP13-FK506 complex is displaced 
by about 2 A RMS relative to the FKBP12-FK506 complex when 
these are overlapped, a point in support of the Yang et al. (1993) 
thesis. 

The structural data presented here for the R42K and R42I mu- 
tant complexes show no such significant rearrangement of the 
40s loop, the 80s loop, or any other part of the FKBP 12 pro- 
tein (Fig. 3). Nor do we observe a change in the conformation 
of FK506, even though we have demonstrated elsewhere (Itoh 
et al., 1995) that just such a conformational transformation is 
present in the FKBP 12 R42K-H87V double mutant complex 
(Kinemage 2). All of our mutant complex structures (including the 
Itoh et at. [19951 double-mutant complex) have been solved in 
the same native FKBP12-FK506 complex crystal form described 
by van Duyne et al. (1991a), a crystal form that includes a sig- 
nificant number of ligand-ligand interactions (van Duync et al., 
1993; Wilson et al., 1995) that might otherwise have compro- 
mised the comparative interpretation weVe presented for these 
structures. Even in the R42K mutant complex, where a signifi- 
cant increase in the temperature factors of the 40s loop of the 
mutant complex structure is observed (see Fig. 4), the weakened 
electron density in this region (Fig. 3C) is still consistent with the 
conformation of the 40s loop found in the native complex struc- 
ture. In the less conservative R421 mutant complex structure, 
however, the fortuitous and unexpected ordering of two water 
molecules (Fig. 3A.B) acting as surrogates of the missing gua- 
nidino nitrogens of R42, leads to a more highly conserved struc- 
ture, even though the loss of CN-inhibitory activity is also 
greater (- 180- fold at 970 nM for the R42I mutant complex ver- 
sus -1 10- fold at 590 nM for the R42K mutant complex; Table 1). 

Elsewhere, we describe the structures of FKBP 12 (Itoh et al., 
1995) and FKBP13 (Griffith JP, Wilson KP, Futer O, Living- 
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ston DJ, Navia MA. Structure of a mutant FKBP13-FK506 com- 
plex that is a high affinity inhibitor of calcineurin [manuscript 
in preparation)) mutant complexes that are inconsistent with the 
hypothesis that significant structural rearrangements in the 80s 
loop region of FKBP12 are responsible (Yang et ah, 1993; 
Clardy, 1993, 1995) for the loss of CN-inhibitory activity seen 
in the corresponding FKBP12 mutant complexes. Nor do the re- 
sults presented here support a similar hypothesis for the 40s loop 
region. Given our structural observations, one might speculate 
that R42 participates in some direct interaction with CN that 
cannot be effectively mimicked by the guanidino nitrogen sur- 
rogate water molecules that stabilize the 40s loop in the R421 mu- 
tant complex structure. With the R42K mutant complex, one 
might further consider an intermediate level of interaction with 
CN, through the single €-amino nitrogen of the lysine substitu- 
ent, and in spile of the greater disorder in the 40s loop. Teague 
(1995) has recently restated the importance of a conformation- 
ally well-defined recognition surface in promoting a strong in- 
teraction between exposed hydrophobic elements, such as are 
found in the FKBP12-FK5G6 effector surface and are presumed 
to exist on the complementary surface of CN. Our results show 
that relatively minor, localized perturbations of the CN com- 
plementarity of the FKBPI2-FK506 effector surface can lead 
to quite significant effects on the CN-inhibitory potential of the 
resulting mutant complexes. 

Methods 

Mutant FKBP12 protein was prepared as reported (Aldape 
et al., 1992; Park et al., 1992; Wilson et ah, 1995). Native and 
mutant FKBP12-FK506 complexes were prepared and crystals 
were grown essentially as described (van Duyne et al., 1991a; 
Wilson et al., 1995). Native crystals were used to seed the mu- 
tant complex crystallization experiments, and all the species re- 
ported here crystallized isomorphously in the native crystal 
form, as shown in Table 2. Diffraction data were collected on 
an XI 000 multiwire area detector (Siemens Analytical Instru- 
ments, Madison, Wisconsin) or on an R-Axis II image plate 
detector (Rigaku/MSC, Woodlands, Texas), as indicated in Ta- 
ble 2. Data collection and processing used software provided by 
the manufacturers. All data were collected at room temperature. 
The reported structure of native FKBP12 in complex with FK506 
(van Duyne et al., 1991a; Brookhaven Protein Data Bank (Bern- 
stein et al., 1977] entry 1FKF) was used directly as an initial 
model for the crystallographic refinement of the mutant and na- 
tive protein complexes. Refinement was by simulated annealing 
using the X-PLOR program package (Brunger, 1992). Mutated 
amino acids were initially refined as alanine, and the actual mu- 
tant side chains were introduced as refinement progressed. Water 
molecules were positioned in the model with the aid of a peak 
search program (SERC Daresbury Laboratory, 1979). The pro- 
gram QUANTA (Molecular Simulations; Burlington, Massa- 
chusetts) was used to examine electron density maps and protein 
models and for the superposition of structures and the calcula- 
tion of the RMS differences reported (Table 2). Coordinates for 
the structures are being deposited in the Protein Data Bank. 
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Note added in proof 

A structure determination of the native FKBP12-FK506 com- 
plex bound to calcineurin has now been reported (Griffith JP, 
Kim JL, Kim EE, Sintchak MD, Thomson JA, Fitzgibbon MJ, 
Fleming MA, Caron PR, Hsiao K, Navia MA. 1995. X-ray 
structure of calcineurin inhibited by the immunophilin- 
immunosuppressant FKBP12-FK506 complex. Cell 52:507-522). 
A preliminary fit of mutant FKBP 12 complex structures to the 
native FKBP 12 in the calcineurin complex is entirely consistent 
with the results presented in this paper. 
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To investigate the ability of a protein to accommo- 
date potentially destabilizing amino acid substitutions, 
and also to investigate the steric requirements for ca- 
talysis, proline was substituted at different sites within 
the long a-helix that connects the amino-terminal and 
carboxy 1- terminal domains of T4 lysozyme. Of the four 
substitutions attempted, three yielded folded, func- 
tional proteins. The catalytic activities of these three 
mutant proteins (Q69P, D72P, and A74P) were 60- 
90% that of wild-type. Their melting temperatures 
were 7-12 °C less than that of wild-type at pH 6.5. 
Mutant D72P formed crystals isomorphous with wild- 
type allowing the structure to be determined at high 
resolution. In the crystal structure of wild-type lyso- 
zyme the interdomain a-helix has an overall bend angle 
of 8.5°. In the mutant structure the introduction of the 
proline causes this bend angle to increase to 14° and 
also causes a corresponding rotation of 5.5° of car- 
boxyl-terminal domain relative to the amino-terminal 
one. Except for the immediate location of the proline 
substitution there is very little change in the geometry 
of the interdomain a-helix. The results support the 
view that protein structures are adaptable and can 
compensate for potentially destabilizing amino acid 
substitutions. The results also suggest that the precise 
shape of the active site cleft of T4 lysozyme is not 
critical for catalysis. 



Phage T4 lysozyme is a small monomeric protein with its 
structure divided into two distinct domains (Fig. 1). The active 
site is located at the junction of the two domains, and it might 
be expected that the alignment of one domain relative to the 
other would be critical for catalytic activity (c/. Storm and 
Koshland, 1970, but see also Jenks, 1969; Knowles, 1991). On 
the other hand, the crystal structure of a fully active mutant 
of T4 lysozyme has recently been described in which there is 
substantial variability in the " hinge-bending angle" between 
one domain and the other (Faber and Matthews, 1990). 

To investigate the ability of T4 lysozyme to compensate for 
disruptive changes in its structure and to determine the need 

* This work was supported in part by National Institutes of Health 
Grant GM21967 and the Lucille P. Markey Charitable Trust, The 
costs of publication of this article were defrayed in part by the 
payment of page charges. This article must therefore be hereby 
marked "advertisement" in accordance with 18 U.S.C. Section 1734 
solely to indicate this fact. 

t Present address: European Molecular Biology Laboratory, Post- 
fach 10.2209, Meyerhofstrasse 1, W-6900 Heidelberg, Germany. 

§ Present address: Laboratory of Molecular Biology, NIDDKD, 
National Institutes of Health, Bldg. 2, Rm. 316, Bethesda, MD 20892. 

% To whom correspondence should be addressed. 



to conserve the alignment of the active site cleft, a series of 
proline substitutions was made in the long interdomain a- 
helix. Studies of proline substitutions and proline replace- 
ments in proteins and peptides include Matthews et al (1987), 
Alber et al. (1988), O'Neil and DeGrado (1990), Strehlow et 
al (1991), and Consler et al (1991), among others. 

The amino acids chosen for substitution with proline, Gln- 
69, Val-71, Asp-72 and Ala-74, are located in the middle of 
the long a-helix (residues 60-80) that connects the two do- 
mains of T4 lysozyme (Fig. 1). It was expected that the 
substitution of a proline at any of these sites would tend to 
significantly distort the a-helix and therefore change the 
alignment of the "upper" and "lower" domain. By making a 
series of replacements it was anticipated that the active-site 
cleft would be distorted in different ways. Also by including 
different substitutions it was possible to include sites that 
were both buried and solvent-exposed. 

Gln-69 is largely exposed to solvent (Table I, Fig. 1) and its 
side chain does not obviously participate in stabilizing inter- 
actions with other parts of the protein. Val-71 is almost but 
not entirely buried. Its side chain makes many contacts with 
the non-polar residues Ile-3, Phe-4, Leu-7, Phe-67, Ala-74, 
He- 100, and Phe-104 that contribute to the hydrophobic core 
both within the carboxyl-terminal domain of T4 lysozyme 
and connecting one domain with the other. Asp-72 is very 
solvent exposed and its side chain is located on the "back- 
side" of the interdomain a-helix relative to the rest of the 
lysozyme molecule (Fig. 1). Ala-74 is largely, but not entirely 
buried, and makes non-polar contacts with residues Val-71, 
Ile-100, Val-103, and Phe-104 that contribute to the hydro- 
phobic core in the carboxyl-terminal lobe of the molecule. In 
wild-type lysozyme Asp- 70 makes an unusually strong salt 
bridge with His-31 (Anderson et al, 1990). Asp-70 also accepts 
a hydrogen bond from the backbone amide of Ala-74. Since 
the replacement of Ala-74 with proline eliminates any possi- 
bility of hydrogen bonding to the amide it was likely that this 
replacement would perturb the Asp-70. . .His-31 salt bridge. 

EXPERIMENTAL PROCEDURES 

Mutagenesis — The proline mutations were introduced by site-di- 
rected mutagenesis according to the uracil template method developed 
by Kunkel (Kunkel, 1985; Kunkel et aL f 1987). The cysteine-free 
"pseudo wild-type" lysozyme (WT*),' in which the 2 cysteine residues 
present in wild-type had been replaced in order to facilitate thermo- 
dynamic measurements (Wetzel et al, 1988; Matsumura and Mat- 
thews, 1989; Pjura et al., 1990) was used as the reference protein. The 
lysozyme e-gene contained within a 630 base-pair BamHl-Hindlll 
fragment had been previously cloned into phage M13 mpl8 yielding 
the derivative M13 mpl8 T4e C54T C97A. It was then transformed 
into Escherichia coli strain CJ236 (dut, ung-, thil, relA/pCJl05(CM')) 

1 The abbreviation used is: WT*, pseudo wild-type. 
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Pig. 1. Backbone of T4 lysozyme showing the locations of 
the proline substitutions discussed in the text. Mutations Q69P, 
D72P, and A74P give active, folded protein. 



which was used to prepare uracil-containing single-stranded template 
DNA. After annealing of the mutagenic primer to WT* template 
DNA, circularization wbb accomplished by using the Klenow fragment 
of E. coli DNA polymerase I and T4 DNA ligase. The double-stranded 
DNA was subsequently used to transform competent E. coli JM101 
cells. Sequencing of the whole lysozyme gene was carried out to 
confirm that no changes had occurred other than the ones introduced. 

Protein Purification— Following subcloning into the plasmid 
expression vector pHN1403 (Muchmore et ai t 1989; Poteete et ai, 
1991) 100 ml of cells grown overnight in 100 ml of LBH broth (10 g 
of tryptone, 6 g of NaCl, 1 ml of 1 M NaOH/liter) was added to 3 
liters of LB broth (12 g of tryptone, 5 g of yeast extract, 10 g of NaCl, 
1 g of glucose/liter) and grown at 37 *C under constant agitation (700 
rpm) and air flow (12 litera/min) in a 5-liter fermenter until the 
optical density at a wavelength of 595 nm reached a value of 1.2. The 
temperature was then decreased to 26 °C, agitation and aeration were 
reduced to 200 rpm and 7 liters/min, respectively. Isopropyl-0-thio- 
galactoside (800 mg) was added to the growth media in order to induce 
lysozyme expression which was allowed to proceed for 100 min under 
continued stirring and aeration. The celts were then harvested into a 
5-liter Erlenmeyer flask where they lysed. A few grains of DNase I 
were added to the now thick and viscous lysate which was left stirring 
at 4 *C for 2.5 h. The lysed cell suspension regained almost the 
viscosity of LB broth, was then taken out to room temperature, and 
placed on a magnetic stirrer for 30 min to allow for complete cell 
lysis. This step almost doubled the final yield of praline-containing 
mutant lysozymes. The lysate became more viscous again and was 
placed back at 4 'C, stirring for another 1.5 h until it regained the 
fluidity of LB broth. All the subsequent purification steps were carried 
out at 4 *C. After centrifugation at 10,000 rpm (17,700 X g) for 2 h, 
only the supernatant contained mutant lysozyme. It was dialyzed 
against 20 mM phosphate buffer at pH 6.5 until the conductivity 
reached a value of 3.8 mS/cm (1 Siemens *» 1 fl The dialyzed 
supernatant was loaded onto a 2.5 X 10-cm CM-Sepharose ion- 
exchange column which had previously been equilibrated with 400 ml 
of 60 mM Tris buffer at pH 7.3. Mutant lysozyme was eluted using 
an 800-ml salt gradient in the range from 50 to 300 mM NaCl in 50 
mM Txib buffer, pH 7.3. Elution was monitored at a wavelength of A 
« 280 nm. Peak fractions (absorbance at a wavelength of 280 nm 
above 0.4 units) were pooled, dialyzed against 50 mM phosphate 
buffer at pH 5.8 for 12 h, and concentrated on a 1 X 5-cm SP- 
Sephadex column previously equilibrated with the same buffer. The 
protein was eluted with SP buffer (100 mM NaP0 4 , pH 6.5, 550 mM 
NaCl, 0.02% NaN 3 ) and stored in this buffer at 4 *C Based on sodium 
dodecyl sulfate and high performance liquid chromatography analysis 



the purity of the mutant T4 lysozyme was estimated to be over 95%. 
Final amounts of 60-120 mg of mutant T4 lysozyme proteins were 
obtained. Wild-type T4 lysozyme typically yields about 150 mg/3- 
liter preparation. 

The activity of the mutant lysozyme was measured at room tem- 
perature using the turbidity assay described by Tsugita (Tsugita et 
al, 1968). Since absolute rates were not reproducible, activities were 
normalized to a wild-type control. 

Measurement of Thermal Stability— Thermal denaturations at pH 
2.0 and 6.5 were monitored by circular dichroism (CD) at a wavelength 
of X = 229 nm on a Jasco J-500C spectropolarimeter as a function of 
temperature (Elwell and Schellman, 1975; Dao-pin et ai, 1990). The 
temperature was varied from 0 to 75 'C at a rate of 1 *C/min with a 
Hewlett-Packard 891 00A temperature controller interfaced to a Hew- 
lett-Packard 87 XM computer. The protein concentration was ad- 
justed to 0.02 mg/ml by measuring the absorbance at \ 280 nm 
using a double beam Varian 2290 spectrophotometer. A probe im- 
mersed in the sample solution just above the UV beam recorded the 
temperature. The solution was continuously stirred with a magnetic 
stirbar. Buffer solutions (150 mM KC1, 10 mM HC1 at pH 2.0 and 150 
mM KC1, 10 mM potassium phosphate at pH 6.5) were prepared from 
doubly deionized, degassed H 2 0 and were filtered before use through 
a 22-jim Millipore filter unit. 

Thermal denaturations were repeated at least three times. Meas- 
urements of the mutants were flanked by WT* thermal denaturations 
under exactly the same conditions in order to minimize errors (Dao- 
pin et ai, 1990). The data were analyzed using standard van't Hoff 
techniques (Becktel and Schellman, 1987; Dao-pin et ai, 1990). 

Crystaliographic Methods — Crystal growth was attempted using 
both hanging-drop 2 as well as batch methods (Weaver and Matthews, 
1987; Alber and MatthewB, 1987) under conditions similar to those 
used for wild-type lysozyme. Both methods yielded crystals of D72P 
isomorphous with WT in the space group P3 2 21. 

Refinement was carried out using the TNT package of refinement 
programs (Tronrud et ai, 1987). The positional coordinates and the 
temperature factors were refined simultaneously using the "conjugate 
directions" option in TNT which improves convergence. 3 

The starting model for refinement was the refined structure of the 
cysteine-free wild-type (WT*) (Pjura et al. t 1990; Bell et ai, 1991) 2 
with residue 72 truncated to Ala. The general approach was to begin 
with low resolution (8-4 A) rigid body refinement with the molecule 
considered as a single unit. This was followed by further rigid body 
refinement but with the mutant molecule divided into two parts 
(residues 1-80 and 81-162). Finally, the molecule was divided into 
three blocks (residues 1-59, 60-80, and 81-162). After examining the 
model on the graphics terminal (Jones, 1978), proline was built at 
position 72 and water molecules from the WT* structure were in- 
cluded in the model. Several cycles of positional refinement using 
moderately weighted geometric constraints were performed using data 
between 20 and 1.9 A. Again, the model was inspected on the graphics 
terminal, some water molecules were added, others repositioned, and 
some side chains adjusted to better fit the electron density. Only 
those water molecules which had a final refined temperature factor 
of less than 80 A 2 , and, in addition, formed hydrogen bonds to protein 
or other bound water molecules and had no steric clashes were 
retained. Thereafter, several cycles of simultaneous positional and 
temperature factor refinement were alternated with model building 
until the crystaliographic residual converged. The number of solvent 
molecules in the refined model is roughly the same as for WT* 
lysozyme, Le. about 145. The refined coordinates have been deposited 
in the Brookhaven Data Bank. 

RESULTS 

Expression of Mutant Proteins 

Gln-69 — ► Pro— Mutant Q69P could be expressed and pu- 
rified in a straightforward manner, yielding up to 100 mg of 
protein from 4 liters of culture medium. The activity of the 
protein is close to that of wild-type (Table I), but it is less 
stable, with melting temperature 12.9 °C less than wild-type 
at pH 2.2 and 7.6 'C lower at pH 6.5. This behavior corre- 
sponds to that of a typical "temperature sensitive" mutant of 

2 A. E. Eriksson, W. A. Baase, and B. W, Matthews, manuscript in 
preparation. 

3 D. E. Tronrud, submitted for publication. 
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Table I 



Sites of proline replacement 



Amino 
acid 


Fraction of side chain 
accessible to solvent 


Mutant 


Relative 
activity 








% 


Gln-69 


0.75 


Q69P 


88 


Val-71 


0.06 


V71P 


a 


Asp -72 


0.76 


D72P 


67 


Ala-74 


0.04 


A74P 


66 



B Purified protein was not obtained for V71P. 



Table II 

Thermal stability of proline-containing mutants 
T m is the melting temperature and &T m the difference between the 
melting temperature of the mutant and that of the pseudo wild-type 
lysozyme (see text). AAG, the difference between the free energy of 
unfolding of the mutant and pseudo wild-type lysozyme, was esti- 
mated from the relationship AAG = AS* AT m (Becktel and Schellman, 
1987) where AS is the entropy of unfolding of the wild- type protein 
(257 and 378 cal/degree mol at pH 2.0 and 6.5). This relationship 
may not be reliable when AT m is large. The quoted values of AAG are 
subject to error, estimated as ±0.5 kcal/mol for Q69P and D72P and 
±1 kcal/mol for A74P. The estimated error in T m is ±0.3 °C for Q69P 
and D72P and ±0.5 °C for A74P. 



Protein 


PH 


T m of mutant 


T m ofWT* 












°C 




kcal/mol 


Q69P 


2.0 


25.6 


38.5 


-12.9 


-3.3 




6.5 


55.6 


63.2 


-7.6 


-2.9 


D72P 


2.0 


28.3 


38.5 


-10.2 


-2.6 




6.5 


56.3 


63.4 


-7.1 


-2.7 


A74P 


3.5 


39.3 


57.5 


-18.2 


-5.7 




4.0 


46.9 


62.2 


-15.3 


-5.0 




5.5 


53.3 


65.7 


-12.4 


-4.5 




6.5 


50.8 


62.9 


-12.1 


-4.6 



T4 lysozyme selected by the random screen of Streisinger et 
al (1961) (c/. Griitter et a/., 1979, 1987; Hawkes et a/., 1984). 
When stored at 4 "C at about 50 mg/ml, Q69P tended to form 
white opalescent aggregates which dissolved on warming to 
room temperature. This process was reversible and seemed to 
have no effect on stability or activity. Small crystals of the 
protein were obtained, apparently non-isomorphous with 
wild-type. 

Val-71 -» Pro— When V71P DNA was transformed into E. 
coli and the ability of a bacterial extract tested to form a halo 
on isopropyl-l-thio-0-D-galactopyranoside lysis indicator 
plates, no halo could be seen at 37 °C. After 24 h at 4 °C, a 
small halo was visible but we cannot rule out the possibility 
that this might be due to a small amount of WT* present as 
an impurity. Attempts to purify V71P by the method de- 
scribed above, or by using a French press to break open the 
bacterial cell walls, or by an alternative method described by 
Dao-pin et al. (1991b) were all unsuccessful. We presume that 
V71P is very unstable and/or is rapidly degraded by proteol- 
ysis. 

Asp-72 — *- Pro— As with Q69P this protein could be readily 
purified by the standard procedure, yielding up to 130 mg of 
protein from a 4-liter culture. Its activity is about 60% that 
of wild-type (Table I) and it is less stable (Table II), again, 
roughly comparable to a typical temperature-sensitive mu- 
tant. Crystals isomorphous with wild-type could be grown at 
4 *C using both batch and hanging-drop techniques. The latter 
method gave the largest crystals, 0.5 X 0.65 X 0.3 mm, from 
2.0 M phosphate at pH 7.1. Some apparently non-isomorphous 
crystals were also obtained at 15 °C but have not been exam- 
ined. 

Ala-74 — Pro— As for Q69P and D72P, A74P behaved 



normally and yielded -140 mg of protein/preparation from a 
4-liter culture. The apparent activity is about two-thirds that 
of wild-type (Table 2). This mutant is less stable than both 
D72P and Q69P (Table II). Because of the low stability the 
protein tends to be partially unfolded at low pH, so the 
stability measurements under these conditions are less relia- 
ble than at higher pH. Some crystals of this protein were 
obtained from 50 mM phosphate, pH 7.9, 16% PEG 8000 and 
are non-isomorphous with wild- type. 

Structure of D72P 

The crystals of D72P appeared to be somewhat more sen- 
sitive to radiation than wild-type lysozyme and did not diffract 
as well. For this reason the exposure time per frame on a 
Xuong-Hamlin area detector system operating with graphite- 
monochromated CuK<, radiation from a Rigaku generator (40 
kV, 100 mA) was increased to 60 s, compared to the usual 30 
s for WT. Under these conditions the data set to 1.9 A was 
89% complete (Table III). 

The map showing the difference in density between D72P 
and WT* lysozyme is shown in Fig. 2a. The positive density 
feature confirms the addition of the pyrrolidine ring. Also 
there is negative density at the site previously occupied by 
the carboxylate of Asp -72. Positive and negative densities also 
indicate the movement of the carbonyl oxygen of Asn-68 away 
from the helix axis (Fig. 26). A negative feature next to the 
carbonyl group of Asp-70 indicates a movement of the oxygen 
toward the helix axis. The side chain of Asp-70 also moves 
slightly, as does His-31 (not shown) maintaining the strong 
(Anderson et aL t 1990) His-31. . .Asp-70 salt bridge. The 
movement of His-31 together with slight adjustments occur- 
ring throughout the lower lobe give the impression that the 
lower part of the long a-helix and the lower domain move as 
an essentially connected unit. 

There is, however, a movement of the upper domain relative 
to the lower one. This shift is shown in Figs. 3 and 4. If 
backbone atoms in the amino-terminal domain of the mutant 
structure (residues 13-59) are superimposed on the corre- 
sponding atoms in WT* lysozyme (Fig. 4a), they have root- 
mean -square discrepancy of 0.16 A, which is essentially ex- 
perimental error. Similarly, the respective backbone struc- 
tures within the carboxyl -terminal domains are also well 
conserved (root-mean-square discrepancy of 0.18 A for resi- 
dues 81-162, Fig. 46). This shows that the structures within 

Table III 
Data collection and refinement statistics 



Data for WT* taken from Eriksson et al} 


Protein 


WT* 
(C54T/C97A) 


D72P 


Data collection statistics 






Mode of data acquisition 


Film 


Area detector 






(Xuong-Hamlin) 


Cell dimensions 






a, b (A) 
c(A) 


60.9 


60.8 


96.8 


98.6 


Resolution (A) 


1.75 


1.9 


Unique reflections 


14,562 


15,147 


Completeness of data (.%) 


65.0 


89.4 


tfror*. (on intensities) {%) 


4.6 


5.5 


Refinement statistics 






Resolution limits (A) 


6.0-1.75 


20.0-1.9 


RMS° deviation from ideal values 






Bond length (A) 


0.015 


0.013 


Bond angle (*) 


2.1 


2.0 


Crystallographic residual {%) 


14.8 


16.7 



° Root-mean -square. 
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Fig. 2. a, map showing the difference in density between D72P lysozyme and wild-type. Amplitudes (F M ut - 
Fwr) where F Mut and Fwt- are the structure amplitudes observed for the mutant and pseudo wild-type crystals. 
Phases calculated from the refined pseudo wild- type structure. Positive contours {solid) and negative contours 
(broken) drawn, respectively, at +3a and -3<r where a is the root-mean-square value of the density throughout the 
unit cell, b, superposition of the A72P mutant structure {solid bonds) on wild-type {open bonds) in the vicinity of 
the mutation. The superposition of the coordinates shown were optimized by least squares prior to drawing the 
figure. 




Fio. 3. a, backbone of the D72P mutant molecule (dark bonds) superimposed on WT* (open bonds). The figure 
is based on the optimal superposition of the amino- terminal region (residues 60-66) of the interdomain a- helix, o, 
superposition of D72P (solid bonds) on WT* [open bonds), as in Fig. 3a but rotated 90V 



the respective NH 2 -terminal and COOH-terminal domains 
are conserved. It is the alignment of one domain relative to 
the other that changes in the mutant structure relative to 
WT*. As shown in Fig. 4, a and 6. this movement corresponds 
to atom shifts up to about 1.5 A. The pronounced maxima 
and minima in Fig. 4, a and 6, are due to different distances 
of the corresponding atoms from the axis of rotation. 

One way to analyze for movements of one part of the 
structure relative to another is to align the mutant structure 
with wild-type based on the superposition of relatively short 
segments of backbone. Two such alignments, based on 6- 
residue backbone segments on either side of the mutation 
site, are shown in Fig. 5, o and b. In Fig. 5a, in which the 
alignment is based on residues 60-66, the amino-terminal 
domain of the mutant structure agrees well with that of WT*. 
This shows that residues 60-66 and the NH 2 -terminal domain 
are connected essentially as a rigid body. In contrast, when 
the superposition is based on residues 74-80 (Fig. 56), neither 



the NH 2 -terminal domains nor the COOH-terminal domains 
of the mutant and wild-type structures coincide. This suggests 
that the residues 74-80 are not connected rigidly to the 
COOH-terminal domain. The upper part of the 60-80 a-helix 
appears to be at least in part responsible for the observed 
flexibility. Small changes in this region can have large effects 
on the rest of the structure. The regions at the beginning and 
at the end of the long helix have been previously identified as 
"hinge-bending" regions in the M6I variant of T4 lysozyme 
(Faber and Matthews, 1990). The hinge-bending that was 
observed for the M6I mutant was assumed to be a low energy 
displacement because different molecules within the same 
crystal displayed different hinge-bending angles. In the case 
of M6I, the long interdomain helix appeared to remain rigid. 
In contrast, in the present case the long interdomain helix is 
bent relative to wild- type. Therefore the hinge-bending seen 
for D72P (Fig. 3) has a very different origin, and presumably 
very different energetics, relative to that seen for M6I. 
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Residue number 

PiO. 4. "Shift plots" showing the displacement between cor- 
responding backbone atoms in the mutant D72P and wild- 
type lysozyme. For each residue the value plotted is the root- mean- 
square discrepancy between the four backbone atoms in the mutant 
(N, C, CA, 0) and the corresponding atoms in WT*. a, superposition 
of the structures of D72P and WT* based on residues 13-59 in the 
amino -terminal domain. The arrowheads indicate those residues that 
are presumed to contact an extended substrate (Glu-11, Asp-20, Thr- 
21, Leu-32, Phe-104, Gln-105, Thr-109, Thr-142 and Arg-145) (An- 
derson et ait 1981). The contacts involving Thr-109 are to sugars in 
subsites A and B, whereas cleavage is between subsites D and E. 6, 
superposition of D72P and WT* based on the carboxyl -terminal 
domain, residues 81-162. 
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FiC. 5. a, superposition of D72P on WT* based on the alignment 
of the amino-terminal part of the interdomain a-helix, residues 60- 
66. 6, superposition of D72P on WT* based on the alignment of the 
carboxyl-terminal part of the interdomain a-helix, residues 74-80. 
Note that a is similar to Fig. 4a indicating that residues 60-66 and 
the amino-terminal domain do not move very much relative to each 
other upon the proline replacement. 6, however, is not very similar to 
Fig. 46, showing that the carboxyl-terminal part of the interdomain 
a-helix does move somewhat relative to the carboxyl-terminal domain 
of the mutant protein. 



Another way of detecting differences between two struc- 
tures is by calculating a difference-distance matrix (Nishi- 
kawa et aL, 1972). All the intramolecular Cf-Cf distances are 
determined for one protein and then compared to the respec- 
tive intramolecular distances of the second molecule. An 
example of such a plot is shown in Fig. 6. Since Cf-C? 
distances in D72P were subtracted from the corresponding 
distances in WT*, positive contours (solid lines) reflect 
shorter distances in the mutant structure; negative contours 
(dotted lines) represent longer ones. As can be seen, the 
mutated residue 72 shows the largest shift with respect to 
residues 95, 122, and 157. 

Fig. 2b illustrates the effect of the Asp -72 — ► Pro mutation 
on the long a -helix itself. In this figure residues 60-66 have 
been superimposed so as to make the structural change more 
obvious. The helix containing the proline has an overall bend 
of about 14 •. A similar analysis of the structure of wild-type 
lysozyme, however, reveals that the same helix already has a 
bend of about 8.5°. Therefore the effect of the Asp-72 -+ Pro 
replacement is to increase the bending by about 5.5°. 

Fig. 7 compares the hydrogen-bond distances within the 
interdomain helix in wild- type lysozyme and in the mutant 
structure. The first a- helical hydrogen bond is from the 
carbonyl oxygen of Thr-59 to the backbone amide of Ala-63. 
Not surprisingly, the introduction of the pyrrolidine ring at 
residue 72 results in a substantial increase in the distance 
between the nitrogen of residue 72 and the carbonyl oxygen 
of Asn-68. The nitrogen -oxygen distances for the successive 
residues are lengthened somewhat, but not excessively so 
(maximum 3.3 A), suggesting that the hydrogen bonds within 
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FiG. 6. Difference distance matrix comparing distances be- 
tween all pairs of a-carbon atoms in the mutant structure with 
the corresponding distances in WT*. The quantity plotted is A,, 
= rij.wT* — rtf,D72P where r^wr* is the distance between the ith and yth 
a-carbon atoms in the structure of wild-type lysozyme, and ry,D72p is 
the distance between the ith and jth a-carbon atoms in the D72P 
mutant structure. Solid contours are drawn at 0.3 A, 0.6 A, 0.9 A r . . . 
and indicate pairs of a-carbon atoms that are closer together in the 
mutant structure than in wild-type. Broken contours, drawn at —0.3 
A, —0.6 A, —0.9 A, . . . indicate pairs of a-carbon atoms that are 
further apart in the mutant than in wild-type. Featureless regions 
indicate domains within which the structure in the mutant is essen- 
tially identical with that in wild-type. 
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FiG. 7. Hydrogen bond lengths within the interdomain he- 
lix. The upper plot shows the distance between the nitrogen atom of 
residue i and the carbonyl oxygen of residue i-4. Within the a- helix 
(i = 63-81) this distance corresponds to a helical hydrogen bond. 
Values for the mutant, D72P, are indicated by the solid line; values 
for WT* are indicated by the broken line. The lower plot shows the 
difference in the distance between the mutant and wild-type struc- 
tures. 

Table IV 



Ramachandran angles for residues within the long a-helix 
of T4 lysozyme 





D72P 


WT* 


A* 




4> 








Ile-58 


-121.3 


169.3 


-126.2 


167.2 


4.8 


2.1 


Thr-59 


-94.3 


171.6 


-90.9 


166.4 


-3.4 


5.2 


Lys-60 


-62.9 


-32.6 


-56.8 


-43.3 


-6.1 


10.7 


Aap-61 


-68.0 


-35.9 


-58.0 


-45.7 


-10.0 


9.8 


Glu-62 


-69.1 


-36.4 


-62.3 


-43.9 


-6.8 


7.5 


Ala-63 


-64.1 


-41.9 


-59.3 


-46.1 


-4.8 


4.2 


Glu-64 


-70.4 


-31.5 


-69.5 


-30.0 


-1.0 


-1.5 


Lys-65 


-61.1 


-49.1 


-69.1 


-42.3 


-8.1 


-6.8 


Leu-66 


-61.3 


-43.2 


-59.4 


-42.0 


-2.0 


-1.2 


Phe-67 


-55.5 


-44.8 


-61.6 


-47.5 


-6.1 


-2.7 


Asn-68 


-62.4 


-22.8 


-56.2 


-43.4 


-6.2 


20.6 


Gln-69 


-89.9 


-38.9 


-64.8 


-40.1 


-25.1 


1.3 


Asp-70 


-66.6 


-38.8 


-67.8 


-37.5 


1.2 


-1.4 


Val-71 


-63.8 


-54.5 


-68.1 


-44.4 


4.3 


-10.1 


Pro-72 


-59.4 


-34.8 


-55.9 


-48.8 


-3.5 


14.0 


Ala-73 


-64.1 


-36.0 


-60.7 


-43.3 


-3.4 


7.3 


Ala-74 


-59.9 


-50.2 


-57.1 


-56.2 


-2.9 


6.0 


Val-75 


-56.6 


-48.0 


-53.5 


-49.7 


-3.1 


1.7 


Arg-76 


-58.0 


-4.20 


-64.4 


-35.2 


6.4 


-6.8 


Gly-77 


-62.5 


-36.3 


-63.4 


-45.1 


0.9 


8.9 


Ile-78 


-64.6 


-46.0 


-57.6 


-46.3 


-7.0 


0.3 


Leu-79 


-65.1 


-17.3 


-69.5 


-20.3 


4.4 


2.9 


Arg-80 


-97.3 


-11.6 


-95.9 


-6.4 


-1.4 


-5.2 


Asn-81 


-93.5 


119.8 


-92.7 


127.3 


-0.9 


-7.5 



the remainder of the a-helix are maintained. In the mutant 
structure a solvent molecule is observed 3.9 A from the 
carbonyl oxygen of Asn-68, suggesting a weak hydrogen- 
bonding interaction. 

The (<t> t yp) angles of the residues within the long helix are 
listed in Table IV. The values observed for the proline (<f> = 
—59.6% $ = -34.6°) agree very well with the average value (<f> 
= —61*, ^ = — 35 p ) for prolines in other protein structures 
(MacArthur and Thornton, 1991). At the site of the substi- 
tution, the addition of the pyrrolidine ring necessitates vir- 
tually no change in <£. The bigger change (AyV = 14*) is in the 
successive peptide. Not surprisingly, the largest changes in 
(0, yp) are for the peptide between residues 68 and 69, which 
is in the previous turn of the a-helix and for which the 
hydrogen bond to the amide of residue 72 is disrupted. 

DISCUSSION 

The most striking result of the present study is the finding 
that proline residues can be substituted at several positions 



within the long interdomain a-helix of T4 lysozyme with only 
modest effects on catalytic activity. The proteins are desta- 
bilized relative to wild-type, but still fold and behave essen- 
tially normally. Attempts were made to substitute prolines at 
four sites and in three cases a functional protein was obtained. 
It is not as if there is one particular site at which a proline 
, can be accepted. Rather, the data suggest that it may be 
possible to introduce prolines at additional sites within the 
interdomain helix, if not at many other sites in the protein as 
well. 

The decrease in stability observed for the two mutants 
D72P and Q69P is very comparable with that found for 
tempera ture-sensitive mutants of T4 lysozyme such as R96H, 
T157I, and A98V identified by the random genetic screen of 
Streisinger et aL, 1961; Griitter et aL, 1969, 1987; Weaver et 
a/., 1989; Dao-pin et aL, 1991a). The mutant A74P is, however, 
less stable than any of these previously described variants. 
The decrease in stability associated with the proline substi- 
tutions seems to be associated to some extent with the inac- 
cessibility of the residue to solvent but the correlation is not 
perfect. Val-71 is largely buried and a proline replacement at 
this site did not yield a functional protein. Ala-74 is also 
largely solvent inaccessible. In addition the proline substitu- 
tion disrupts the a-helical hydrogen bond between residues 
70 and 74 which, in turn, is likely to misalign and perhaps 
weaken the very strong salt bridge between Asp-70 and His- 
31 (Anderson et a/., 1990). In this case protein was obtained 
although with substantially reduced stability. O'Neil and 
DeGrado (1990) found that the energy cost of an alanine to 
proline replacement within a dimeric a-helical model peptide 
was 3.4 kcal/mol. Also Yun et aL (1991) found, by free energy 
simulations, exactly the same value for an alanine to proline 
replacement within a short polyalanine helix. These values 
are roughly comparable with those found here (average values 
of 3.2, 2.7, and 5.2 kcal/mol) but there is no reason to expect 
close agreement since the context is different in every case. 

In the case of A72P, for which the crystal structure is 
available, the pyrrolidine ring is seen to introduce three 
unfavorable contacts with neighboring atoms (2.93 A between 
C 4 and the peptide nitrogen of Val-71; 3.04 A between C 4 and 
the peptide nitrogen of Val-71; 3.00 A between C 4 and the 
carbonyl oxygen of Asn-68). Each of these contacts corre- 
sponds to an unfavorable van der Waals interaction energy 
in the range 1-2 kcal/mol (Levitt, 1974). These steric clashes 
are, presumably, a major factor in the destabilization of the 
mutant structure relative to wild-type. 

The results provide further evidence that protein structures 
are adaptable and can compensate for amino acid substitu- 
tions at many sites (Matthews, 1987; Sondek and Shortle, 
1990). It also illustrates the redundancy that is present in the 
amino acid sequence of a protein. Not every amino acid in 
the linear sequence is necessary for folding (Reidhaar-Olson 
and Sauer, 1988; Zhang et aL, 1991). Amino acid substitutions 
that are expected to distort and destabilize the folded struc- 
ture of the protein, and disrupt a major a-helix that might be 
a key folding intermediate, do not prevent the formation of a 
folded functional protein. 

Functional proteins were obtained with prolines substituted 
at positions 69, 72, and 74. Residues 69 and 72 are located in 
successive turns on the same side of the helix. Residue 74, 
however, is on the opposite side of the a-helix. Therefore one 
cannot argue that prolines are only accommodated on one 
side of the a-helix such that each substitution bends or 
distorts the a-helix in the same direction. 

The substitution of a proline in an a-helix that is part of a 
protein is not the same as a proline substitution of an isolated 
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a-helix. In the absence of any constraints, a substituted 
proline is likely to completely disrupt the helix (Strehlow et 
al t 1991) or to introduce a bend of up to 45° (Barlow and 
Thornton, 1988; Karle et oi., 1991). The structural conse- 
quences of a proline introduced into an a-helix in a protein, 
however, will be modulated by interactions between the a- 
helix and the rest of the protein. Prolines in a-helices in 
known protein structures are associated with bend angles that 
average 26 ± 5° (Barlow and Thornton, 1988), but this is not 
to say that an introduced proline will cause a change of this 
magnitude. The determination of the structure of D72P pro- 
vides an example of the structural compromise that can result* 
The presence of the proline increases the bend angle of the 
a-helix, but only by 5.5*. Thus, it is not to be expected that 
each individual proline substitution will cause a major re- 
arrangement of the two domains in T4 lysozyme. Nevertheless 
it can be anticipated that the different substitutions will cause 
distinct changes in the alignment of one domain relative to 
the other. This suggests that the precise shape of the active- 
site cleft in the resting enzyme is not critical for catalysis. At 
the same time it should be noted that the key catalytic residues 
are at the base of the active site cleft, and the structural 
changes in this region may therefore be relatively small (<0.5 
A). The residues in T4 lysozyme that are presumed to contact 
an extended oligosaccharide substrate are indicated in Fig. 
4a. 

It should be noted that the crystals that were obtained for 
the replacement Asp- 72 — > Pro are isomorphous with wild- 
type. It is possible that the formation of such isomorphous 
crystals constrains the structure of D72P to be more similar 
to WT* than is the case in solution. The differences between 
the structures of D72P and WT* in solution may therefore be 
larger than those seen in Figs. 3 and 4. The mutants Q69P 
and A74P give crystal forms that are not isomorphous with 
wild-type. Hopefully, the structure analysis of these crystals 
will give a better overall impression of the structural changes 
that are induced by the proline substitutions. 

The first determination of the structure of a temperature- 
sensitive mutant lysozyme (Grutter et a/., 1979) showed that 
relatively large changes in stability were accompanied by 
minimal changes in structure, these being localized to the 
immediate vicinity of the substituted amino acid The ability 
of T4 lysozyme to accommodate destabilizing mutations in 
this manner has subsequently been seen on a number of 
occasions (Matthews, 1987) and is further exemplified by the 
present study. 
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